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ABSTRACT 
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A new model for spatial distributions of reflectivity in storm cells has been fitted to 
digital radar data from the National Severe Storms Laboratory (NSSL) in Oklahoma for 
May 1968 and from the MIT Lincoln Laboratory in Massachusetts for July 1970. The 
NSSL data were taken with a modified WSR-57 weather radar with 2.6-km resolution. 

The NSSL data consisted of modified B-scan records on magnetic tape of storm cells 
tracked at 0° elevation for several hours. The Millstone Hill MIT L-band radar with 0.8-km 
resolution produced cross-section data on several cells at 1/2° elevation intervals. The model 
developed uses ellipses for contours of constant effective-reflectivity factor Z with constant 
orientation and eccentricity within a horizontal cell cross section at a given time and eleva- 
tion. The centers of the ellipses are assumed to be uniformly spaced on a straight line, with 
areas linearly related to logZ. This asymmetrical model can display the steep gradient in 
log Z at some cell aspect found by Shaw and others. All cross sections are similar at different 
heights (except for cell tops, bottoms, and splitting cells), especially for the highest reflec- 
tivities; wind shear causes some translation and rotation between levels. Goodness-of-fit 
measures and parameters of interest for 204 ellipses are considered. The model was found 
to be a good representation of the reflectivity and attenuation of intense storm cells, although 
more investigation is required of the problem of cell splitting. 
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Section 1 
INTRODUCTION 


1.1 GENERAL 

This document is the final report on the storm modeling study performed by Com- 
puter Sciences Corp. (CSC) for NASA under contract NAS5-21638. The general problem is 
prediction of interference between space and terrestrial communication systems caused by 
precipitation scattering. Specifics of the problem and the approach to its solution are sum- 
marized in the following section, which includes a summary of previous work by CSC in the 
field and an indication of work to be done. Other sections present radar data, methodology, 
storm modeling, results, and conclusions. It may be noted that the storm models developed 
have application to studies of attenuation of radio waves as well as scattering, and also to 
various meteorological storm studies and hydrological simulations of the water cycle. 

1.2 SUMMARY 

To support the sharing of frequency bands by space and terrestrial radio communica- 
tion services, the CCIR (International Radio Consultative Committee) has instituted Study 
Programme 5- ID- 1/5, including the determination of the average distribution in time and 
space of observed values of effective scattering cross section per unit volume of precipitation 
as a function of height above ground, time of day, season, and climatic region. The observa- 
tions used in the present study are digitized reflectivity-factor (Z) data from radars at the 
Massachusetts Institute of Technology Lincoln Laboratory (MIT-LL) and the National Se- 
vere Storms Laboratory (NSSL). The data are essentially horizontal sections through storm 
cells. (See table 1-1.) 


Table 1-1 .—Radar and Data Characteristics 


Source 

NSSL 

MIT-LL 

Area 

Oklahoma 

New England 

Date 

May 1968 

July 1970 

Time 

1735 to 2020 

2310 to 0003 

3-dB beamwidth 

2.0° 

O 

o 

Range resolution 

2.6 km 

0.8 km 

Tilt 

0° 

1° to 5° 


1 




The approach used was to assume a model, fit it to the data, and compute the error of 
fit. The basic model for a contour of constant reflectivity is an ellipse of the same area, as 
suggested by Rogers and Rao (ref. 1) as well as Amorocho (ref. 2). The former authors as- 
sumed ellipses all of the same orientation, with those in a set for a given cell and height con- 
centric but with different eccentricities. The model used here assumes different orientations 
for different sets, but the same orientation and eccentricity within a set. Furthermore, the 
ellipse centers of a set are assumed uniformly spaced on a straight line. As suggested by pre- 
vious work and that of Hudlow (ref. 3), the square roots of the areas within a set are con- 
sidered a linear function of log Z. This asymmetrical model can display the discontinuity in 
reflectivity at some cell aspect found by Shaw (ref. 4) and Donaldson. 1 

The assumed model was found to represent the contours of a storm cell cross section 
at a given height very well in general. A median error E of 25 percent was found in fitting 
individual ellipses to 204 contours: 


E = 


100 


U- X 
U 


where 

U = union, area within either ellipse or contour 

X = intersection, area within both ellipse and contour 

The median additional error resulting from the constraints of the above assumptions on 
ellipses within a set was less than 5 percent. It was found that three horizontal sets, for 
heights of about 5, 6, and 7 km, could be well represented by the same set model, so that 
the proposed cell model consists of nested, vertical elliptical cylinders, displaced or rotated 
at some heights by wind effects, with rounded or plumed tops and flatfish bottoms. A typi- 
cal horizontal cross section, with logZ considered as a third dimension, is an asymmetrical 
elliptical paraboloid. 

In spite of displacement and rotation of the basic cylinder at some heights, the peak 
reflectivity Z Q and the slope of the fitted linear profile, log Z 0 /( max area), varied very 
slowly spatially from height to height with the MIT-LL data, or from time to time at a given 
height with the NSSL data. The samples were too small to establish appropriate storm cell 
classes or characteristic parameters for each. However, it appears that the data resolution 
was adequate and that the model is accurate and general enough to describe the intense cells, 
which seem to be the most definitely structured. The chief problem encountered was incip- 
ient or actual splitting of a single cell at one or more levels. 

A comparison was made of the attenuation along paths through the centers of five hori- 
zontal sections of one cell (MIT-LL cell A) as observed and as modeled by ellipse sets. The 
agreement at three of the levels was excellent, but at two levels it was degraded because the 
orientation of the core log Z = 4 ellipse remained relatively constant with height, whereas 
the average rotated counterclockwise (ccw). It may therefore improve the model for some 

1 Ralph J. Donaldson (Air Force Cambridge Research Labs.), personal communication, Nov. 1970. 
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purposes to use the core orientation for the ellipse set instead of the average of the individual 
orientations. 

Many more storm cells should be processed to verify the model parameters derived and 
to suggest appropriate methods for treating splitting cells, or at least to indicate the magni- 
tude of the problem. 
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Section 2 

OVERVIEW OF STORM MODELING 


! Preceding page 


2.1 INTRODUCTION 

In section 1 the broad problem was stated to be prediction of precipitation-scatter 
interference. More specifically, it is to predict distributions of effective scattering cross sec- 
tion per unit volume as functions of time and space in two scales with different orders of 
magnitude. Permissible interference is defined in terms of very small fractions of an hour, 
and precipitation is known to vary considerably in distances less than a kilometer. So on the 
one hand, distributions with resolutions of seconds and meters are of interest; on the other, 
these distributions are wanted as functions of season and climate, with scales of months and 
hundreds of kilometers. It is not surprising that the solution of this problem is not yet at 
hand. 

2.2 APPROACH 

One possible type of solution is experimental, the measurement of enough known con- 
figurations in enough places for sufficient time. As will be shown, this is impractical even 
in a rather trivial case. A more fruitful approach is analytic, to study the nature of the dis- 
tributions of interest with a view to relating them to some broad existing data base. The 
only available data useful for this approach are radar reflectivity and surface rainfall. Only 
weather radars specially equipped for digital storage can provide useful resolution in time 
and space, including height, and useful quantities of storm records. (See sec. 4.1.1.) Thus, 
for a first approach, Rayleigh scattering is assumed and radar backscatter data are used. 
Second, the only data base broad enough in the large scales of space and time is found in the 
world’s rainfall archives of approximately a hundred years. Not every country, however, can 
match the U.S. hourly precipitation data from gages spaced about 50 km (30 n. mi.) apart in 
many areas. Fortunately, the United States is subject to a wide variety of climates, so that 
results found should not be difficult to extrapolate to similar climatic regions elsewhere. 
Therefore, the problem is resolved into a study of the characteristics of reflectivities aloft 
using radar data, of rainfall on the surface, and determining their correlations so that one 
may be predicted from the other. 

Studies of reflectivity and rainfall require reduction of large volumes of data to a small 
number of characterizing parameters. This process can be based on assumption and evalua- 
tion of models. (See sec. 4.1 .) These models may ultimately be somewhat sophisticated, like 
the one for storm cells proposed and given preliminary evaluation in this report. A simple, 
yet powerful model developed in earlier work is the linear storm profile stressed in this sec- 
tion. In general, a profile is a plot of one variable as a function of another, but such a curve 
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may require many parameters or the plot itself for storage and manipulation. It is much more 
convenient if appropriate scales can be found for the quantities of interest so that the profile 
is, with acceptably small error, a straight line. Two numbers, a slope and intercept or two 
intercepts, then characterize the function for efficient storage and simplified usage. The 
finding of such useful models is stressed. It will be noted by reference to the literature that 
only models involving extensive averaging in time and space, and hence of limited applica- 
tion, have been available, and that those found in the previous work (secs. 2.4 and 2.5) and 
only recently confirmed by others are in some ways still ahead of the state of the art. Two 
applications seem to merit particular attention. The first is the use of a field of rain gages 
as a sort of low-resolution radar with hourly scans to provide area sampling that, from pre- 
liminary results, would avoid the pitfalls of sampling sparse data at one point only and pro- 
vide profiles analogous to those found with radar data. The second interesting application 
of the linear profile is to use radar data quantized to a few levels to characterize about a 
60-dB range in reflectivity, to derive peak values with a useful degree of accuracy. 

2.3 EXPERIMENT VERSUS ANALYSIS 

It would seem that storm models could be developed either by experiment or by study, 
or perhaps by an optimal combination of studies checked by experiments. However, there 
is an inherent danger in relying on experiments based on insufficient sampling of rare events. 
This is well brought out by the reluctance of a British expert, E. G. Bilham, to rely on even 
50 years of point-rainfall records. He indicates that even decades without intense rains in a 
small area are more likely due to chance than to climate and states (ref. 5) that “the really 
intense rains lasting no more than an hour or two seem to be distributed without regard to 
the annual average rainfall.” The principle involved is not, of course, to avoid experiment, 
but rather to avoid reliance on too little data to warrant the conclusions desired. It will be 
seen below that this problem is overcome by use of radar or rainfall data from a circle of 
about 185-km (100 n. mi.) radius; i.e., area rather than point sampling. 

2.4 SUMMARY OF PREVIOUS WORK 
2.4.1 Composite Radar Profiles 

The earliest extensive body of ra : nfall reflectivity data was taken at X-band with an 
AN/CPS-9 radar at Montreal Airport by McGill University’s Stormy Weather Group from 
May to September 1963. CSC Report 453340 (ref. 6), prepared under NASA contract 
NASW-1216, used this data base and appropriate geometry, often approximate, to develop 
the preliminary results for intersecting beams that are reported in CCIR Report 339. 

Monthly profiles included all types of rainfall and the underlying assumption was made that 
the reflectivity was uniformly distributed; i.e., equally likely to occur anywhere in the 185-km 
(100-n. mi.) circle of observation or in the worst month of observation. In spite of some ob- 
vious possibilities for error, this was a useful first approach. Unfortunately, McGill’s present 
radar is not appropriately equipped to produce digitized data, and early radar data of NSSL 
were not satisfactory, although some of those data are considered below and more should 
be processed for additional information. 
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2.4.2 Hourly Radar Profiles 

Even though the early NSSL data were not accurate, the errors were apparently sys- 
tematic so that many relationships were preserved; the data were used in CSC Report 6012-1 
(ref. 7) to derive characteristics of height and reflectivity profiles which have proven useful 
as bases for later work. The data used were taken by NSSL using a modified WSR-57 radar 
at Norman, Okla., from April to June 1966. As detailed in reference 7, the S-band data were 
quantized in square units 2.6 km (1.4 n. mi.) on a side in steps of 10 dBZ (dBZ = 10 log Z) 
and were used to produce profiles of height versus the square root of the fractional area 
covered by effective reflectivity factor Z above some value, and of log Z versus the same 
abscissa. On the scales indicated, the profiles were surprisingly linear and constant in slope 
for several consecutive hours. More or less well-separated high- and low-slope categories 
were found and tentatively labeled for convenience as convective or continuous, respectively. 
Useful correlations were found between the slopes and other parameters of the height and 
log Z profiles and in particular between the intercepts of the log Z profiles and the simul- 
taneous radar tops (highest echoes within about 185 km (100 n. mi.) at a quarter before 
each hour) observed at the_Oklahoma City radar of the National Weather Service network. 
These findings gave hope that a commonly assumed Z-R relationship, where R = rainfall rate 
in millimeters per hour, would provide log Z profiles from surface data, so that with the 
heights derived from radar tops and use of two classes of slopes or modes, improved predic- 
tions would be possible. 

Accordingly, the next stage of the work concentrated on surface rainfall, but two other 
matters should be noted before leaving the radar study. First, correlations were found, in a 
separate study, between radar tops all over the United States and hourly rainfalls having an 
expected 25-yr return period, as found in standard publications. These correlations provide 
bases for useful criteria as to minimum permissible heights for intersections of beams of 
space and terrestrial communication systems. Second, selected storm cells from similar 1 968 
NSSL data were used for development of an ellipse-fitting technique, to be covered in more 
detail. It was found that the linear profiles of the type mentioned were spurious in convec- 
tive cases involving only single or closely similar cells. The simplified early method of line 
fitting had obscured a maximum-height cutoff found by another method and recognized as 
applicable to single cells. 

2.4.3 Rainfall Profiles 

In CSC Report 6012-2 (part of Report RC-10245 to the FCC; ref. 8) several types of 
rainfall profiles were considered, using standard hourly precipitation data from 30 stations 
fairly uniformly distributed within a 220-km (120-n. mi.) radius of Norman, Okla. It was 
shown that straight lines generally resulted from plotting cumulative log R distributions ver- 
sus the square root of fractional time, or area, or area-time product, such as for 30 stations 
during 20 hr in a month. Such area sampling in terms of N, the number of stations reporting 
more precipitation than R mm hr - 1 , exhibited a smoothing tendency explained by two- 
dimensional extension of the Nyquist sampling theorem and the relatively stable nature of 
the storms as they pass over one station after another. Not only were linearity and stability 
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found as with radar data, but also the high and low classes of slope. A study of the literature 
found agreement with the results of Huff (ref. 9) and Hudlow (ref. 3) as to the suitability of 
the log/? to the square root of area relationship for durations up to a few hours. However, 
Huff used a gage spacing of only about 5.5 km (3 n. mi.) instead of the ordinary widely 
spaced data used here. Also, his smaller network with sides 37 km (20 n. mi.) long did not 
completely cover the widespread low-slope continuous rain configurations. In Report 6012-2 
(ref. 8) one storm that evidently consisted of only a single cell in the whole circle of gages 
was found to have the same linear profile against area found from single-cell radar data. Be- 
cause of the poor sampling resolution, this profile was relatively inaccurate. The profile was 
not checked nor compared with the radar data for the same period. 

2.5 STORM CELL CHARACTERISTICS 

The elemental storm unit to be modeled is the cell. A basic postulate was assumed from 
observation that even though rainfall may be found in any and all configurations, the high- 
intensity, high-extent cells of most importance result from organized activity, such as in 
the supercell discussed by Browning (ref. 10) and others. Thus order may be found where it 
is of most interest. The ellipse was assumed for the basic shape of a contour of constant re- 
flectivity, in the horizontal plane, as was done by Rogers and Rao (ref. 1), providing more 
flexibility than the circles used by Dennis (ref. 11). The assumptions used here regarding rela- 
tionships between contours are related to those of Rogers and Rao as shown below: 


Assumption 

Rogers and Rao 

Altman 

Elliptical shape 

Yes 

Yes 

Same eccentricity 

No 

Yes 

Same orientation 

All 

All in a set 

Same centers 

Yes 

No 


Although Rogers and Rao assumed concentric ellipses, in this report the centers of the 
ellipses in a set for different values of Z at a given time and elevation are assumed to be uni- 
formly spaced on a straight line. This is in agreement with the data and results of Shaw 
(ref. 4) and of Amorocho (ref. 2). 

The remaining factor required to define the model is the relationship of the area within 
a contour versus log Z at the contour, vhich gives linear profiles. It was noted previously 
that the general relationship found in CSC Report 6012-1 was linear for log Z versus the 
square root of the area. However, for single cells, the abrupt changes of log Z normally found 
at one aspect of the cell in the raw data, and confirmed by meteorologists experienced in the 
field, could only imply linearity of log Z versus area, or length squared, a parabolic function. 
This was also found, as noted, with rainfall data in the case of a single-cell storm. Therefore, 
it has been concluded that isolated cells or complexes of closely similar cells exhibit profiles 
linear with area, and that aggregates of dissimilar cells exhibit profiles linear with square root 
of area. The synthesis of an aggregate from individual cell profiles has not been, but should 
be, attempted. 
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The data bases used for the present study of cell characteristics have been given in table 
1-1 . Results from NSSL data were obtained by the author in a previously unreported study. 
The study reported here, based on a small amount of data from these two sources, showed 
the utility of the model assumed and validated its correctness for a set of ellipses at a given 
time or tilt angle. But the limited data provide only preliminary notions as to the nature of 
time or space ensembles. It appears, reasonably, that the ellipse parameters vary only slowly 
from time to time at a given elevation, but vary considerably with height in a single cell be- 
cause of the change with height of the winds aloft. At present, the best cell model appears 
to be one of substantially vertical elliptical cylinders, with portions displaced and rotated 
somewhat by wind shear, and with rounded tops and flat bottoms. Configurations at cell 
birth, death, and splitting will require more data and study. 

2.6 WORK TO BE DONE 

Now that radar reflectivity, surface rainfall profiles, and cell characteristics have been 
studied separately, it is important to establish the relationship between peak surface rainfalls 
and peak reflectivities aloft. This relationship— 

(1) Is complex because rain may not fall until late in the storm life cycle and will then 
be strongly localized in time and space 

(2) Has been obscured by the usual averaging in time and space 

(3) Should be studied with profile technique using both slopes and intercepts, which 
is more powerful than using averages alone 

NSSL radar data and Agricultural Research Service (ARS) rainfall data from the same area 
in Oklahoma are at hand for April and May of 1969 and available for the same period of 
1970. The radar outputs have a resolution of 1 .85 km (1 n. mi.), with tilted antenna data to 
indicate vertical structure. The ARS network comprises about 200 stations on about 5.5-km 
(3-n. mi.) centers, providing an unusually large, dense system. The profile techniques already 
developed, and described above, should be used to establish detailed space and time correla- 
tions between peak surface rainfall and peak reflectivity aloft to facilitate prediction of pre- 
cipitation scatter interference. 

As the CCIR sharing criteria permit harmful interference no more than about 0.01 per- 
cent of the worst month, it is clear that rare events such as storms associated with floods are 
included, so that methods used in their prediction should be of value. In fact, Bilham and 
other U.K. and U.S. hydrologists use statistics of extreme values. It is not generally realized 
that extremes can be more valuable predictors than medians. Further study in this area may 
be profitable. 

Several items suggested by work done in previous studies were noted previously and are 
summarized here: 

(1) The linear log R to the square root of the area cell profile found in CSC Report 
6012-2 (ref. 8) should be checked for accuracy against radar data and the exercise 
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repeated for other similar data to establish the utility of surface data from normal 
widespread low-resolution gage networks. If the ARS-NSSL comparison proposed 
above is performed, the high-resolution ARS network can be progressively de- 
graded by omission of alternate stations to produce even more valuable results. 

(2) The synthesis of aggregate profiles from individual cell profiles should be studied 
to confirm the compatibility of the two profile relationships found in previous 
studies. 

(3) The data base of cell characteristics should be extended far beyond the few cells 
considered here, which are too few to support any sound conclusions. Because it 
has been found that useful accuracy can be obtained from coarse data by the pro- 
file technique, several years of NSSL April-May data should be processed to insure 
that the assumptions used here apply to at least the most intense cells and in un- 
usual circumstances. The incidence of splitting should be noted, and if it is found 
statistically frequent, attempts should be made to model it usefully. Cell sorting 
could be automated, and statistical summaries of parameters of most interest 
should be programed. 
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Section 3 

DIGITIZED RADAR DATA 


3.1 INTRODUCTION 


Before describing in detail the two types of digitized rad/ar data used for this study, it 
is necessary to establish a uniform nomenclature. This is done with reference to figure 3-1. 

A header line carries a date-time group and other items, followed by a data sequence, a term 
used loosely here. The data lines are called records , made up of boxes, each containing a 
single-digit value of reflectivity factor Z. As the digit/may represent Z linearly or logarithmi- 
cally, the term “intensity” is used for the digit in general. The contiguous group of records 
making up a cell, or occasionally two cells or a ebnfiguration about to split or otherwise 
irregular, is called a block. A group of blocks for a given time and tilt, such as the top two 
in the figure, is designated a quadrangle ox quad. 


3.2 NSSL DATA 


The weather radar data system at the NSSL has been described by Wilk et al. (ref. 12). 

It is based on the WSR-57, the standard weather radar of the National Weather Service, which 
has the following characteristics— 


Wavelength 
Peak power 
Pulse length 

Pulse repetition frequency 
Minimum detectable signal 
Beam width 
Antenna gain 
Azimuthal scan rate 


10 cm 
450 kW 
4 jus 
164 s' 1 
-110 dBmW 
2 ° 

38.5 dB 

2 rpm for PPI displays 

0.4 rpm for magnetic tape recording 


It has been modified by addition of digital video processing equipment including the Lher- 
mitte integrator developed at NSSL. Lines with reflectivity of log Z = 4 or greater from the 
data recorded on magnetic tape, provided by Kathryn Gray of NSSL, are displayed in fig- 
ure 3-2. Each sequence starts with a month-day-year-time group and initial range, incremen- 
tal range, and tilt angle in degrees. The reflectivity data follow in rotated B-scan format 
(bearing vertical, range horizontal). The first three digits of each line are the azimuth of the 
2°-wide antenna beam. The following 77 digits 0 to 5 represent the average logZ in “boxes” 
2° wide and 2.624 km (1 range unit(RU) or 1.417 n. mi.) long. The 77 RU’s extend from 
26.28 to 228.4 km (14.19 to 123.3 n. mi.). Angles are computed on a grid with spacings of 
1° by 1 RU, so they are true at about 1 50 km (80 n. mi. or 50 RU), near the center of the 
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range. On the printout grid, distances in the azimuth direction are true at about 125 km (67 
n. mi. or 41 RU). It should be noted that the B-scan or polar grid is treated as rectangular, 
as this is closely true for the small cells considered. The raw data are culled to reduce com- 
putational load and noise by removing— 

(1) Fixed echoes within certain azimuths and ranges 

(2) Azimuths containing only 1 ’s 

(3) Strings of l’s and single 2’s, isolated within an azimuth 

(4) Isolated azimuths 

(5) Strings of l’s and single 2’s, isolated within a range slot 

The resulting preprocessed data then appear as in figure 3-1, ready for ellipse fitting, with a 
header for each quad and a separator after each block to define its range limits. 

3.3 MIT DATA 

The MIT data are from the Millstone Hill L-band radar of MIT-LL, furnished by Dr. 
Robert Crane. Table 3-1 tabulates the radar characteristics and table 3-2 shows factors 


Table 3-1.— MIT-LL Radar Characteristics 


Frequency 

1 .295 GHz (23.2-cm wavelength)' 

Antenna 

26-m (84-ft) parabola with Cassegrainian feed 

Antenna gain 

46.7 ± 0.3 dB (“new” subreflector) 

Beamwidth 

0.7° between half-power points 

Polarization 

Right-hand circular transmitted 
Left-hand circular received 

Transmitted power 

3.3-MW peak (continuously monitored) 

Pulse length 

12.4 jus 

Pulse repetition rate 

120 s' 1 

Receiver bandwidth 

80.5 kHz (12.4-ps matched predetection filter) 

Resolution volume 

1° by 1.9 km 

Data processing 

Analog-to-digital conversion of sine and cosine 
channels every 1 0 ps 

Computer sampling rate 

20 s' 1 

Detection 

Square law by computer operations 

Dynamic range 

80 dB accomplished by combining the output 
from two receivers 

System noise temperature 

280 K (includes atmospheric and ground effects 
averaged over 0° to 30° elevation angle) 

Overall system line losses 

2.8 dB for receiver 1 
43.5 dB for receiver 2 

Matched filter processing loss 

1.1 dB 

Single-pulse value for unity signal-to-noise ratio 

-30 dBZ at 100 km 
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Table 3.2.— MIT-LL Measurement Accuracy 


Measurement 

Absolute 

accuracy, 

maximum, 

dB 

Calibration 

repeatability, 

maximum, 

dB 

Sample 

fluctuation, 

rms, 

dB 

Transmitter: 




Power monitor 3 

0.2 

0.1 

0.1 

Line loss 

.1 

— 

— 

Antenna: 




Transmit and receive gain 

.6 

— 

- 

Pattern integration 

.1 

- 

- 

Receiver: 




Noise tube 

.3 

.1 

- 

Calibrate source: 3 




Power monitor 

.2 

.1 

- 

Attenuator setting 

- 

.1 

- 

Line loss 

.1 

- 

- 

Estimation error for Rayleigh process 

- 

- 

.7 

Measurement of received-to-transmitted 




power ratio 

.7 

.2 

.7 


includes line loss error (maximum) used to transfer readings to reference points. 


affecting radar measurement accuracy. Figure 3-3 outlines the region scanned southwest of 
the radar in the hour just before midnight of July 29, 1 970, during a thunderstorm. The 
quads of interest are indicated, as well as typical blocks. Echoes in part of the lower left-hand 
quad and part of the one above it are shown with dBZ values in figures 3-4(a) to (i). The 
headers show radar tilts or elevation angles from 1° to 5° in 1/2° steps at times from 2310 
to 2357 hr. The echoes of most interest are marked A, B, Cl, and C2 on figure 3-4(c), show- 
ing a scan where all appear. 

The computer preprocessing of the MIT-LL magnetic tape record is summarized in fig- 
ure 3-5. Briefly, the original tape is reblocked as it is copied and quad maps of log Z are pre- 
pared. From the latter, blocks of most interest are selected. These are mapped as in figure 
3-6. The boxes are 0.7813 km on a side, with intensities given as log Z to one significant 
figure, so that a 4 is used for 4 < log Z < 5. The line or record numbering is arbitrary. On 
the corresponding punched cards it is multiplied by 2 for compatibility with the ellipse-fit- 
ting program and the 2° azimuths of the NSSL data. Individual cells in a block are selected 
and extraneous points are rejected. As shown in figure 3-5, maps and cards are prepared with 
intensities quantized on linear thresholds as well as logarithmic. The linear thresholds are 
based on the maximum block intensity M as follows: 
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Digit Fraction of maximum intensity M 


5 

4 

3 

2 

1 

0 


0.8 to 1 
.6 to .8 
.4 to .6 
.2 to .4 
.1 to .2 
Oto .1 


Arbitrary values of M can also be inserted when a storm cell appears on two quads. 
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Figure 3*1.— NSSL data blocks. 
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Figure 3-3.-MIT-LL radar scan. 
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-73.66 


DATE 72970 TIME 2322BL0CK 
DELTA X .781 DELTA Y .781 
QUAD ORIGINS, Y< -160.00, 


ORIGINS, Y< 


-73.10 


- 120 . 00 , 


1 0000000000000000000000000000000000000000000000000000 
2 0000000000000000000000000000000000000000000000000000 

3 0000000000000000000000000000000000000000000000000000 

4 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 0000000000000000000 

5 0000000000000000000000000000000000000000000000000000 

6 0000000000000000000000000000000000000000000000000010 

7 oooooooooooooooooooooooooooooooooooooooooooooooooooo 

8 000000000000000000000000000000000000000000001 1100000 
9 0000000000000000000000000000000000000000001111100000 

10 00000000000000000000000000000000000000001 11 110000000 
11 00000000000000000000000000000000000000001 11 100000000 
12 0000000000000000000000000000000000000000010000000001 

13 0000000000000000000000000000000000000000000000000000 

14 0000000000000 00000 0000000000000000000000000000 11 00 11 

15 0000000000000000000000000000000000000000000011000111 

16 0000000000000000000000000000000000000000000000000110 

17 0000000000000000000001000000000000000000000000000000 

18 00000000000000000001 11000000000000000000000000000000 

19 0000000000000000001110000000000000000000000000000000 

20 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOdOOOOOOOllOOO 

21 0000000000000000000000000000000000000000000000111000 
22 0000000000000000000000000000000000000000000000110000 

23 0000000000000000000000000000000000000000000000000000 

24 0000000000000000000000000000000000000000000000000000 

25 0011000000000000000000000000000000000000000000000000 
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Section 4 

STORM MODELING METHODOLOGY 


4.1 INTRODUCTION 


4.1.1 General 

Storms have been characterized by profiles that have been produced and used by hy- 
drologists for many years, as discussed quite generally by Court (ref. 13). These profiles are 
usually averaged over many hours and are not accurate enough for more than rough calcu- 
lations. As Eagleson (ref. 1 4) says on page 1 96 of Dynamic Hydrology: 

Except in the most simple cases, it is usually quite difficult to deduce an ac- 
curate lagrangian description (i.e., relative to the translating fluid) of a mesoscale 
disturbance from an array of eulerian (i.e., fixed-point) measurements. However, 
this lagrangian description is essential to the further understanding of storm struc- 
ture and to the improvement of generalized methods for hydrologic forecasting. 

. . . Weather radar studies of storm structure should ultimately provide the desired 
lagrangian data. 

As noted in the previous section, weather radar data are being used for this study. How- 
ever, it is not just an accurate description that is being deduced but a model— a simplified 
description expressing a large quantity of data with a few parameters preserving the essential 
characteristics of the storm cells. The development of such a model can be considered as 
several separate problems. 


4. 1 .2 Contour Model 


The basic data configuration is the contour around the area within which the reflectiv- 
ity is higher than some level, say log Z- 4. The first problem, which is not considered here, 
is to isolate such contours from other contours of equal interest and from those of lesser 
interest which represent only computational noise. The second is to select the basic shape to 
fit a contour of interest. Dennis (ref. 1 1 ) used the circle as a useful first step. Court sug- 
gested (ref. 13) the ellipses of the bivariate normal distribution, followed recently by 
Amorocho (ref. 2). Rogers and Rao (ref. 1 ) used ellipses, as will be described. A sample of 
1 5 1 ellipses fitted to log Z data is discussed in section 7.3.1 , and the median error is found 
to be near 20 percent. Thus the selection of the ellipse seems justified. It also seems clearly 
preferable to the circle in view of the elongated shapes that were found, with minor/major 
axis ratios as low as 0. 1 . 




* hi h' • v 
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4.1.3 Set Model 

Fitting an ellipse is a two-dimensional problem, but each storm cell cross section or 
sequence may have up to five contours for different values of log Z, so fitting these may be 
considered a three-dimensional problem, with log Z as the third dimension. One objective 
is, of course, to find simple relationships between the ellipses to reduce the total number of 
descriptive parameters. The bivariate normal distribution suggested by Court is a possibility, 
but its spread at low values does not fit the data, so Amorocho (ref. 2) used only the upper 
portion. Rogers and Rao (ref. 1) also assumed concentric ellipses of constant orientation, 
but with variable ellipticity. The model assumed here and justified later uses constant ellip- 
ticity and orientation, but ellipse centers are uniformly displaced along a straight line. Two 
functional relations between Z and area have been tried, log Z versus area and log Z versus 
square root of area. 

4. 1 .4 Four-Dimensional Models 

The above models may add a fourth dimension in two ways considered here. The NSSL 
scans at 0°-tilt angle are repeated at intervals of 1 5 min. MIT-LL scans for several tilt angles 
in a single cell involve log Z or Z as a function of the usual three dimensions. Thus ellipsoids 
could be fitted, but the variations of wind with height would shear the ellipsoids so that a 
piecewise approach is preferable. It is found that several tilt sequences (roughly the same as 
heights) may be considered as identical, forming a family of approximately concentric ellip- 
tical cylinders with a different family above some height. Another reason for not using el- 
lipsoids is the anomalous nature of the top and bottom of a cell, frequently with a plume at the 
top and an irregular bottom. One unsolved problem that should be mentioned is that of splits. 
Cells often split into two after developing a second nucleus, assuming intermediate config- 
urations that are not well treated by arbitrary division into one or two cells. At the least, it 
must be decided by characterizing many more cells whether this phenomenon is statistically 
significant. 

4.2 FITTING ELLIPSES 

The problem of finding the best fit line to an irregular cluster of points, (x p y.), i = 1, 
n, is an old one in the field of statistics. Usually the regression line of y on x, y = ax + b, is 
found. Suppose the regression of x on Y, :: = ay + b', is also computed. By rearrangement 

1 b' 

y= x . 

a a 

leading one to expect a - l /a' and b = ~b'la', which is generally not the case. This is of little 
concern in the usual applications of regression lines, but when ellipse axes are being sought 
that will (in some sense) best fit the scatter of points, more symmetry is required than the 
regression lines are capable of giving. 

It is known that any pair of coordinate axes (u, v ) will exactly split the scatter plot 
into equal halves if their origin is the center of mass. The problem is determining which co- 
ordinate system will equalize the points among its quadrants. It is a fact that such a 
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coordinate system cannot have a regression line as one of its axes, so a different approach 
must be taken. 

Consider the covariance matrix 




0 11 

0 12\ 

v°21 

hJ 

'2>? 

Zv«\ 


u) 

o 12 . For an elongated 


be more points in the first and third quadrants than in the second and fourth. Hence, for 


a l2 =1 x iVi 


there will be more products of like sign (yielding positive terms) than of unlike sign (yield- 
ing negative terms). Therefore, a x 2 =£ 0. Clearly the coordinate system ( u , v), to equalize 
the points in its quadrants, will have ct' 2 = 'Lu ( v j = 0. 

Expression of (w ; ., u ; .) in terms of a rotation of (x /; y.) yields 

u j = x ( . cos 6 + y . sin 6 

v i = ~ x i s * n ^ + y i cos ® 

Substitution into the expression for T gives the rotated matrix 


r'= 


r ' = ( 


°'l2\ 

\ 

V CT 21 

°2 2/ 





Vit- 

!•?/ 

j x cos 2 6 + a 2 2 sin 2 6 + o 12 

sin 2 6 

Vi(o 2 2 “ 0 j j ) 

x h(o 22 - Uj j ) sin 26 + o l2 cos 26 

Uj 1 sin 2 6 + a. 


Since a' 21 = a'j 2 , the 6 value needed to make o' 2l =0, makes o' 12 = 0 also. The degree 
of symmetry required is an orthogonal rotation of axes to equalize quadrants. This approach 
also works in three or more dimensions, because the covariance matrix is always symmetric. 
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Of course, a suitable 6 may be found by solving for a [ 2 = 0, or 


6 = — arctan 



But this technique is useful only in two dimensions. In more dimensions at least three angles 
are involved, and the formulas become much more cumbersome. To gain insight into a more 
general method, matrix algebra may be used. If two matrices A and B are related by 

A = PBP r 

where P is an orthogonal matrix and T denotes transposition, then determinant A = deter- 
minant B; A and B have the same eigenvalues; and the eigenvectors of B and A are identical. 
In particular, if A is diagonal, then the eigenvalues of B are the diagonal entries of A and the 
eigenvectors of B are orthogonal. For the scatter plot considered, the rotation of axes is an 
orthogonal transformation and the ( u , v ) axes lie in the unique positions of the widest and 
narrowest spread of the scatter plot. Hence, the question arises of whether this relationship 
holds between T, T', and the orthogonal rotation matrix 

( cos 6 sin 9\ 

-sin 6 cos 9/ 


The answer is affirmative, because a matrix multiplication will verify that 


r' 


= FtrR r 


Now eigenvector techniques may be brought to bear on the problem. A method that 
will produce the eigenvalues and eigenvectors of a matrix will give the a! . values and the 
actual axes of the coordinate system that equalize the quadrants. The o' n values, being vari- 
ances, indicate the spread of the points along the new axes. To find the axis of the widest 
spread, the eigenvector of the largest eigenvalue (largest o ' H ) is used. To find best-fit ellipses, 
semimajor and semiminor axes are chosen proportional to the values a! (standard 

deviations) and with values such that n times their product equals the area being fitted. 


An illustration of the above remarks is afforded by the bivariate normal distribution 
suggested by Court (ref. 13) and applied by Amorocho (ref. 2). Its density is given by 

(1/2 Tto x o y y/\ - p 2 ) exp j - [( x/o x ) 2 - 2p(x/a x )(y/o y ) + (y/o,) 2 ] } 

_____ 

where the correlation coefficient p = cr /a a . Points of equal probability density lie on 

x y x y 

ellipses 


K = 



■ 2p — — + 

oo 
x y 



— - 2a — — + — 
a x y a o o 
xx xx yy yy 
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where K > 0. These ellipses have the same shape as an ellipse with semiaxes a x and a , and 
their orientation 6 with respect to the (x, y ) coordinate lines may be found from these 
formulas. 


The error of fit for two figures of the same area must consider shape. The measure used 
originally was (A 0 /A^) X 100, where A f is the figure area and A q is the area of the boxes 
with centers outside the ellipse. It was later found that a more general mathematical concept, 
the symmetric difference metric, which is twice the above for equal areas, was used for this 
purpose by Lee and Sallee (ref. 15). Thus, [1 - (Ay- A Q )/(A e + A Q )] X 1 00 is now used for 
generality, where A g is the ellipse area. 

4.3 FITTING SETS OF ELLIPSES 


The individual-ellipse fits to storm cell data as displayed at the left of figure 4-1 indi- 
cate the reasonableness of using the average of all thetas (TAV) and of all ratios of minor/ 
major axes (RAV) for all ellipses of a given set. It is also seen that the centers are not far 
from uniformly spaced on a straight line. This line may be found from a three-dimensional 
version of fitting a line to a group of points so that the sum of the squares of the errors is 
minimized, assuming that log Z or Z is the third dimension. As the points being fitted are 
the centers of gravity or centroids of the individual ellipses, the resulting line is called the 
centroid line. It is specified by the three block coordinates of a point on the line and two 
angles. 


It remains to examine how ellipse area varies with log Z. If log Z is considered as a third 
dimension, the problem may be visualized as finding the shape of the appropriate solid. It 
will first be shown analytically that the two simplest alternatives are the cone and the 
paraboloid. 


Assuming for simplicity that the ellipses have their major axes along the y axis, their 
equations are 



a 2 b 2 

But mb = A = area and ajb = R = axis ratio, so 



: 2 +R 2 v 2 =-=— 


IT 


K l 


where 


K 2 =7T. 


33 



If the square root of the area is a linear function of log Z 

A = (B - S log Z) 2 


where 

B = intercept of line with\/A axis 
S = slope of line 

then the cross section in the plane y - 0 is 

±Kx = B - S log Z 

This is a linear function, so the surface is a cone. Similarly, if the area is a linear function of 
log Z 

a = b - s log Z 

then the cross section at y - 0 is 

K 2 x 2 =b - s logZ 


This is a parabola, so the surface is a paraboloid. 

Both of these alternatives are displayed in figure 4-1 . The paraboloid at the center dis- 
plays two features of interest: First, it does not come to a point like a cone, but has a blunt 
top which may not quite reach the log Z = 5 level. This is not important because the equa- 
tions are available for scatter or attenuation calculations up to the fitted peak of perhaps 
log Z = 4.9. It is probably more accurate to determine the shape of the top from the main 
portion of the data and not to place too much weight on the two boxes with log Z = 5. 
Second, and most important, the paraboloid models the asymmetry found in nature by 
means of its crowded contours along one edge. It may be noted that the two lowest contours 
touch in both the individual fits and in the paraboloid. The preference for this model is rein- 
forced by discussions with meteorologists such as A. J. Chisholm and Ralph Donaldson 
studying severe storms and by data such as appear in figure 4-2, taken from Shaw (ref. 4). 
The outline of the cloud is shown as a thin line; the heavy lines within the shading (denoting 
radar-detected precipitation) are isopleths of Z = 2 X 1 0 2 , 2 X 1 0 4 , and 2 X 1 0 6 mm 6 m _ 3 . 
Outlines are shown at time t of 10, 30, and 50 min after the storm first appears. 

For the results reported in section., 5 and 7, the few points with log Z = 5 were ignored 
to avoid the possible large quantizing errors associated with the smallest areas, and the area 
profiles were offset as required to pass through the points at log Z> 4. This action was in- 
tended to give stronger weight to the higher reflectivities. It was found later, however, as 
reported in section 6, that comparisons of attenuations on paths through the centers of the 
actual cells and the corresponding models indicated that the offset and omission of the 
points for log Z = 5 increased the errors. The result is believed to be critical only for this 
measure and the central paths. It may also be advisable in this case to set the common ellipse 
orientation to that of the highest level instead of to the average of all. For other uses, these 
details may not be significant. 
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4.4 FITTING RELATED SETS OF ELLIPSES 


The sets of ellipses considered so far have each been derived from data taken almost 
simultaneously at a given time and tilt or elevation angle. The question now arises of rela- 
tions between sets derived from the same cell but made at different elevation angles or at the 
same elevation angle but at different times. These are clearly related and storm modeling 
should demonstrate this stability in time and space, or any major lack of continuity such as 
caused by wind shear with height or by splitting. Conceptually, this part of the study is also 
curve- or line-fitting as above, such as noting constancy of orientation, or linearity of motion. 

In line with this discussion, the only effort made to fit related sets of ellipses was for 
different elevations in the same storm cell. It would be of advantage to model the three-di- 
mensional cell as not just a stack of unrelated ellipse sets but as an elliptical cylinder or an 
ellipsoid. Alternatively, the approach may be taken of finding the errors resulting if one of 
these assumptions is made. The computer program was made interactive so that similarities 
could be forced between the ellipse sets at different elevations to show these effects. The sig- 
nificance of the results is difficult to interpret. It was found, for example, that the same 
ellipse sets can be used at elevations from 2° to 4.5°, with differences only in translation 
and rotation of complete, otherwise identical, sets at the different levels. This translation 
and rotation may have a great deal of interest for the meteorologist, but would not be sig- 
nificant for a horizontal radio ray from a random direction. It is perhaps advisable to gather 
more data before formulating more hypotheses or conclusions. 

The interactive procedure for relating sets of ellipses is as follows. The data pertaining 
to all levels of a cell must be isolated, stored, and culled. Using appropriate switches on the 
computer console, a preliminary run may be made with or without ellipse plotting, as de- 
sired. Then on a second round, some levels may be given the same centroid lines as others 
that are closely similar. The same line at different locations and orientations may be tried 
for all the echoes between selected levels. On a last run, the area sets may be equalized at 
all levels. On all runs the parameters including error are tabulated and plots made. Results 
are presented in the next section. 
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Figure 4-1.— Typical ellipse fits to NSSL data. 



Figure 4-2.— McGill reflectivity contours. 
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Section 5 

EXAMPLE OF MODEL FITTING 


5.1 INTRODUCTION 

Because of the complexity of the model-fitting techniques, and because several sources 
and types of data are involved, it seems wise to present the results in two separate ways. In 
this section, log Z data from one storm cell, MIT-LL cell A, are considered in detail to clarify 
the methodology and the terms and formats used. In the next section, other cells, linear Z 
data, and NSSL temporal sequences are considered, and standardized results are compared 
to bring out similarities and differences. 

The basic results are all presented in the form of computer printouts. The abbreviations 
used in the printouts are explained in table 5-1 and a typical printout is shown in table 5-2. 
The printout format is as follows, keyed to table 5-2. 

(1) An overall header contains a date group and labels for the five columns corre- 
sponding to the contours of minimum logZ = 1, etc., as displayed in the table. 

(2) The tabulation of data for each time is headed by the time and data relating to 
the individual contours. 

(3) Below the header extending across a whole line, more data for individual ellipses 
are shown for comparisons with data for related ellipses, in following lines. 

(4) The first data line starting at the left margin describes the ellipse family auto- 
matically fitted using the assumptions of section 2.5. 

(5) Following lines under the same time header note the results of changes of many 
possible types, generally with the object of showing uniformity between succes- 
sive related sets. 

The same time header may be repeated for additional cells from the same block, or for 
splits. In this case, data for related sets become separated. 

5.2 FITTING SETS OF ELLIPSES 

An example of the raw logarithmic data to be fitted has been shown in figure 3-6. It 
provides four contours, around the digits 1 or greater, 2 or greater, etc., as shown at the top 
of figure 5-1. These are fitted with ellipses with the results printed on the top lines for 
time = 2322 of table 5-2 under the headings Min Z = 1, 2, 3, and 4. To help in following 
further processing, these parameters are plotted on the X, Y plane in figure 5-2(a) with the 
semimajor and semiminor axes at the correct orientation of the former, THTA. Note that the 



Table 5-1 .—Abbreviations for Tables 5-2 and 7-3 to 7-5 


Abbreviation 

Definition 

MODAYR 

Date-MOnth DAy YeaR 

MINZ=1 

Refers to contour including log Z > 1 

TIME 

Time 0 to 2400 

THTA 

Theta-azimuth of upper semimajor axis, ccw from north in 
degrees (MIT); subtract AZ for NSSL 

MAJ 

Length of semimajor axis in boxes (MIT) or RU (NSSL) a 

MIN 

Length of semiminor axis in boxes (MIT) or RU (NSSL) 

RAT 

Ratio of MIN/MAJ 

CFA 

Centroid azimuth -90° (NSSL) or Y in boxes (MIT) 

CFR 

Centroid range in RU (NSSL) or X in boxes (MIT) 

CFZ 

Centroid log Z 

TRZ 

True azimuth of centroid line projection (MIT); subtract CFA 
for NSSL 

PRJ 

Projected spacing between log Z = / and log Z = / + 2 

TAV 

Average theta or true azimuth (MIT); subtract CFA for NSSL 

RAV 

Average ratio of MIN/MAJ 

CG-AZ 

Center-of-gravity azimuth (NSSL) or Y in boxes (MIT) 

CG-RU 

Center-of-gravity range (NSSL) or X in boxes (MIT) 

AREA 

Ellipse area in RU deg (NSSL) or square boxes (MIT) 

ER 

Error of fit in percent (see sec. 4.2) 


a Range units (RU) = 2.624 km (1.417 n. mi.); 1 box = 0.7814 km. (See sec. 3.2.) 


scale for the axes is not the same as for the rest of the figure. The centers of gravity, or cen- 
troids, are projected on the X-\og Z plane in figure 5-2(b), showing the dashed centroid line, 
CENLN, fitted for least-square error. PRJ gives a measure of the slope of this line, the pro- 
jection of the CENLN between X, Y planes log Z = I and log Z = I + 2. The CENLN is shown 
projected on an X, Y plane in figure 5-2(c). The nominal c.g. coordinates for Min Z = 1 to 4 
and 2.5 are printed in the first full line for time = 2322, but the values are truncated by drop- 
ping hundredths, so they do not agree precisely with the figure. The CENLN azimuth, TRZ, 
is shown as well as the average THTA, TAV, common to all major axes. The areas within the 
contours, from the same printout line, are plotted versus Min Z in figure 5-2(d). A profile is 
fitted for least-square error, as shown, - T 'ith the apex marked at the maximum log Z, for 
which the area is zero. The slope is arbitrarily defined as 

maxIogZ 

max area 


= 100 


4,75 

140 


100_ 100 _ 
DA 29.5 
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Table 5-2.— MIT-LL Ellipse Data Printout 
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where DA = differential area per unit log Z. The last items printed are the errors of fit. To 
give this unconventional measure discussed in section 4.2 some physical significance, the 
original contours and the fitted ellipses have been superposed in figure 5-1 . The computed 
errors (in percent) are— 


MinZ 

1 

2 

3 

4 

Cell A 

23 

25 

11 

8 

Cell B 

17 

28 

10 

37 


Table 5-3 summarizes not only some of the data already presented, but also the aver- 
ages and standard deviations of the fits for MIT-LL cell A. The figures must be interpreted 
with some care. For RAV, the average ratio of minor to major ellipse axes, a coefficient of 
variation of 0.057/0.70, or about 0.08, is obtained. For TAV, however, the base should not 
be the average azimuth, but 360°, for a coefficient of about 0.05. For the centroid lines, the 
rms error is in linear box units of 0.78 km, and for the area in square box units. The uniform- 
ity of many of the parameters, such as TAV, RAV, and max log Z, is quite remarkable if 
data from erratic top and bottom levels are ignored. 

5.3 FITTING RELATED SETS OF ELLIPSES 

In figure 5-3 the sets of fitted ellipses for horizontal cross sections at all heights of cell 
A are shown. The columns from left to right are raw data contours, individually fitted el- 
lipses, and ellipse sets as described above. The fitting errors are collected in table 5-4. Only 
two of the individual fit errors are over 30 percent. Min Z = 4 at 2345 has note a in the 
table and figure. It is seen that the high error is due to the ragged shape of the contour and 
to its small size, which leads to high quantizing error. Note b shows a high error from a 
split at the cell top. Set fit errors more than 1 0 percent higher than the individual fit errors 
are also noted. Notes c, e, and g refer to Min Z = 4 cores maintaining the same orientation 
at all cell heights even though the lower reflectivity structures rotate increasingly ccw with 
height. Notes d and f seem due to a random orientation of Min Z = 3 contours. 

Figure 5-3 illustrates the effects of changes forced on the ellipse parameters by the 
interactive procedure. The method is in general to average parameter values from all levels, 
and to force all levels to the same value. This reduces the number of parameters required 
to describe the cell and leads to the simplified elliptical cylinder model. The first, least 
significant, change was to set all CENLN projections PRJ equal to 0.70 as indicated at the 
head of a column of the figure and at one line in the table by the abbreviation P = .70. As 
expected, the ellipse centers are now all uniformly spaced, and the only significant error 
increase is at note h because of the large original projection. The common value of 
MIN/MAJ = RAV = R = .65 caused no error increases over 10 percent. Forcing a common 
major axis orientation 65° right of north produced large error increases at notes j to n, only 
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Table 5-3.-MIT-LL Ellipse Set Data 
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Table 5-4.-MIT-LL Ellipse Fitting Errors 
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k occurring at Min Z = 4. In the right-hand column of the figure, R = .65, T = -65°, and 
the area profile has been chosen for best fit for time 2322 to 2340 as 

Area = 180- 40 log Z 

This caused large errors at notes p, q, and r at the top and bottom of the cell, as expected, 
but it is seen that the proposed model of sheared, nested, elliptical cylinders is validated. 



Figure 5-1.-MIT-LL fitted ellipses. 




SCALE FOR MAJOR 
AND MINOR AXES, 
RANGE UNITS 

I 1 

0 5 





Figure 5-2.— MIT-LL plotted ellipse data, (a) Individual ellipses, (b) Centroid line fit. 
(c) Ellipse set fit. (d) Area profiles. 
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Section 6 

TEST OF STORM MODEL 


6.1 TEST METHOD 

For the development of the storm model, error of fit was estimated using geometric 
criteria as discussed above. Its intended application, however, involves transmission of elec- 
tromagnetic energy, so a more relevant criterion of usefulness is desirable. Therefore, the 
attenuation of a radio wave through the observed reflectivity pattern has been compared 
with that through a storm cell reconstructed from the model parameters. This is facilitated 
by using east-west or north-south paths parallel to the X and Y axes of the MIT-LL data. 

For these paths, precise values of reflectivity are available for each increment of 0.78 km. 

It will be shown for the three east-west central paths through five cross sections of the 
MIT-LL cell A, that Z sums by the two methods generally differ by about 5 percent. 

6.2 CALCULATION OF ATTENUATION 

With certain usual assumptions about the dropsize spectrum, the attenuation A in deci- 
bels through a storm cell may be expressed as 


A = K Jz a ds 


where 

K= a function of frequency, dB km -1 
a = constant to be discussed 
Z = reflectivity, mm 6 m -3 
ds = distance along path, here 0.78 km 

For greatest accuracy, the Constanta should be considered to vary as in table 6-1, taken 
from Zawadski and Rogers (ref. 1 6), but if the values of K shown are used, a may be taken 
as unity without large errors, for Z values of most interest, around 1 0 4 . The last two rows 
show the attenuation corresponding to a typical Z ds of 1 0 5 and the attenuation error for 
5 percent error in the integral as a function of frequency. 

6.3 RESULTS 

A computer program was written to construct digitized horizontal cross sections of 
cells similar to those in figure 3-1 from the fitted ellipse parameters, and to sum the values 
of Z in the rows in the X direction. When these sums were compared with similar sums 
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Table 6-1.— Attenuation at Various Frequencies 


Attenuation factors 

Frequency, GHz 

10 

20 

30 

40 

Best a 

0.78 

0.70 

0.66 

0.63 

Best K for a= 1 

2 X 10“ 5 

10" 4 

3 X 10 -4 

5 X 10 -4 

A for / Z £&=10 5 , dB 

2 

10 

30 

50 

0.05 A, dB 

0.1 

0.5 

1.5 

2.5 


Table 6-2.— Highest Values of f Z ds/ 1000 


Time 

From data 

From model 

Error for 

highest peak, percent 

2322 

149 

199 

125 

181 

206 

129 

3 

2328 

87 

129 

139 

88 

122 

132 

5 

2334 

219 

348 

266 

168 

217 

203 

38 

2340 

83 

111 

72 

63 

70 

63 

37 

2345 

43 

50 

42 

52 

53 

43 

6 


derived from the raw data, it was found that agreement was improved by using area profiles 
without the offset to produce equal areas at the log Z > 4 level, which is described in sec- 
tion 4.3. It was also found essential to use the point at log Z = 5 in fitting the area profile 
for time 2334. In table 6-2, the paths with the three highest sums are compared for MIT-LL 
cell A at all elevations at which log Z > 4 was observed. The errors are high at 2334 and 
2340 because the ellipses for log Z > 4 are not alined with the others in the sets. If north- 
south paths were considered as well as east-west, and the highest integrals of both sets con- 
sidered, this error would be reduced. In any case, most of the errors are small at low 
frequencies. 
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Section 7 
RESULTS 


7.1 INTRODUCTION 

The storm modeling results presented are based on data of three types: 

( 1 ) MIT-LL elevation sequence— logarithmic 

(2) MIT-LL elevation sequence— linear 

(3) NSSL time sequence— logarithmic 

The first type provides spatial characteristics in three dimensions from successive radar 
sweeps at increasing elevation angles. The second considers distributions of reflectivity Z 
instead of log Z to focus on the high-intensity cell cores. The third type provides temporal 
characteristics by following several cells for several hours in two dimensions with a horizon- 
tal radar scan. 

To facilitate comparisons, all three types are considered in the same order used previ- 
ously; data first, then fits to individual ellipses, sets of ellipses, and finally related sets of 
ellipses. 

7.2 DATA 

7.2.1 Logarithmic MIT-LL Data 

Before considering results, the data used will be reviewed and displayed in the more 
convenient form of figure 7-1 . Here the portions of figure 3-4 containing cells A and B have 
been used in the form of digital maps like figure 3-6. Contours have been drawn about the 
l’s, 4’s, and 5’s, etc., and placed side by side for comparison. There is a common Y grid, the 
displaced X references are labeled, and the horizontal scale is marked. Cell A, in .the upper 
half of the figure, displays intensity levels 1 to 4 from 2° to 4° in elevation, or from about 
5 to 9 km in height. Only a single log Z = 5 box appears, at 3°. The bottom tapers down sym- 
metrically but the top shows a split at the log Z = 2 level. The centers of the 1 0 to 3‘/2° levels 
are close to a vertical line; higher levels have drifted to the east. The time required for scan- 
ning is not negligible, so the 36 min intervening between the 1° and 4° scans may be 
significant. 

Cell B, also shown logarithmically, does not display much vertical uniformity. In fact, 
a split is indicated by the peninsula at the bottom at 2° and is definite at 1 V 2 0 . Only the 
northern part shows at 1°. More and faster scans such as those of NSSL for 1969 and later 
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are required to separate fully the effects of time and elevation in a case like this. Without an 
extended time sequence, the two effects may easily be confused. 

Another aspect of splitting is evident in cell C, shown logarithmically in figure 7-1. 
Definite splits exist at both the top and bottom. Small internal peaks at IV 2 0 to 2V° are dis- 
regarded, and the 1 Z° and 2° levels are split near the peninsulas, and consistent with this at 
2 V 2 0 as shown. The part called C2 is fairly weak and uniform from W 2 0 to 3°. The other 
part; Cl, probably continues below 1°, or a height of about 3 km. 

7.2.2 Linear MIT-LL Data 

The cell A linear data are plotted below the logarithmic cell A data in figure 7-1. The 
significance of the shades used in the figure may be seen from table 7-1. The right-hand 
column applies for time 2334. In the logarithmic data, the digit 4 is used for 40 < 10 logZ< 50, 
or in this case, for 40 to 50, whereas the linear presentation uses five levels to provide greater 
detail in the small intense area responsible for most scattering and attenuation. 

7.2.3 NSSL Data 

NSSL data are used as well as MIT data because the former include long time sequences 
giving preliminary information on cell velocities, lives, splits, etc. The automatically blocked 
data of figure 3-1 were manually culled to remove trivial echoes, and significant related 
echoes were identified as shown and grouped for best display of their characteristics in fig- 
ure 7-2. Cell scans each 15 min 'oetween 1735 and 2020 of May 3, 1968, were studied. 
Elongated cell B split into cells B and C between 1820 and 1835. 

7.3 INDIVIDUAL ELLIPSE FITS 
7.3.1 Fitting Errors 

The results of the individual ellipse fitting as described above are included in tables 7-3 
to 7-5 (presented subsequently), which use abbreviations explained in table 5-1. To facilitate 
understanding the results, these are extracted and presented in various ways. The first ques- 
tion of interest is the accuracy of the fits. Using the error measure discussed at the end of 
section 4.2, figure 7-3(a) characterizes the fit errors in a histogram format giving the number 
of cases in each error category. The emrs over 50 percent are due to splits, as are many of 
those from 40 to 50 percent. 

Data sets for two NSSL storm cells are shown in figure 7-2 with the left-hand sets of 
ellipses showing the individual fits. For example, in the upper cell there are two 5’s, one 
above the other, so there is a vertically oriented central ellipse. The 4’s ellipse is sharply tilted 
to the right, and the 3’s and 2’s ellipses have a slight tilt. The l’s have not been used because 
they are insignificant electromagnetically, pruned arbitrarily, and often in a very ragged con- 
figuration, as in a plume at the top of a storm. The apparent misalinement of the 5’s ellipse 
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Table 7-1.— Linear Scale Equivalents 


Shade 
(fig. 7-1) 

z/z 

' max 
linear 

101og(Z/Z max )X 10 

Typical dBZ 

5 

0.8 to 1 .0 

9 to 10 

50 

4 

.6 to .8 

8 to 9 

49 

3 

.4 to .6 

6 to 8 

47,48 

2 

.2 to .4 

3 to 6 

44,45,46 

1 

.1 to .2 

0 to 3 

41,42,43 

0 

0 to .1 

0 

<40 


ccw of the average azimuth and of the 4’s clockwise (cw) is probably due to quantizing 
error. That is, there are too few samples (boxes too large— low resolution) to define the con- 
tour accurately. 


7.3.2 Cell Sizes 

As a preliminary estimate of cell size for the ellipses fitted to the contours of log Z = 4 
for the MIT logarithmic data and the contours Z = 0.6Z max for the MIT linear data, the aver- 
age semiaxis lengths are as follows: 


Cell 

Semimajor 
axis length 

Semiminor 
axis length 

MIT A (logarithmic) 

2.9 

1.8 

MIT B 

2.7 

1.8 

MIT C 

4.5 

2.5 

MIT A (linear) 

2.0 

1.1 

MIT C 

2.2 

1.7 


As the lengths are given in units of 0.78 km, a typical cell at the log Z = 4 level has major 
and minor axes of about 4.5 and 2.8 km, respectively. The geometrical mean, or radius of 
the equivalent circle, is 3.5 km, as noted in various CCIR reports. The cell sizes can be found 
for other levels, but the highest levels tend to be the most irregular, thus it is well to derive 
such data from the lines fitted to several levels, or sets of ellipses, as will be discussed. The 
many other individual cell characteristics may properly be ignored after it is shown that they 
are adequately represented by the set parameters. 
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7.4 SETS OF ELLIPSES 


7.4.1 Fitting Errors 

Figure 7-4 graphically portrays the individual ellipse parameters for MIT-LL cell A as 
in figure 5-2(a). Inspection of the diagram for the different elevation angles shows both regu- 
larity and irregularity. The diagrams for 2.5°, 3°, and 3.5° are very similar, taking account 
of the expanded scale. Those for 1 .5° and 2° appear somewhat similar until it is noted that 
the 2° centroid line points in the direction opposite to those of the other four. The 4° and 
4.5° diagrams form a small rough group, but it seems reasonable that the bottoms and tops 
of the cells would be most irregular. To emphasize this, the log Z = 4 ellipses at 2.5°, 3°, 
and 3.5° are closely identical, whereas the 1° and 5° levels shown in figures 3-4(a) and 3-4(i) 
are highly irregular. Fortunately, they are also very weak, an inverse of the basic rule postu- 
lated previously that intensity and order go together. Similar plots for the fitted sets of el- 
lipses are shown in figure 7-5, but only one ellipse axis is shown for each set, as they all have 
the same orientation. 

The errors of fitting the ellipses to the basic contours are shown later in tables 7-3 to 
7-5. It is of interest to study the increase in error between the individual fits and the set fits. 
These are shown in histogram form in figure 7-3(b). The differential error due to set fitting, 
that is, alining orientations and eccentricities, and fitting centroid and area lines has a me- 
dian value of about 5 percent. This remarkably low figure indicates the basic orderly struc- 
ture of the cell, even though the individual contours are not precisely elliptical. It is strongly 
implied that the basic assumption of elliptical shapes is reasonable. 

7.4.2 Set Characteristics 

The characteristics of the ellipse sets for a given time or elevation are best considered 
along with those for the related sets, but a few points are worth noting separately. For ex- 
ample, the cell size at one-half peak reflectivity may be found by averaging the areas from 
table 7-4 for Min Z- 3, as this corresponds to Z/Z max = 0.5 according to table 7-1. Attention 
is restricted at first to those times when a peak of 40 dBZ occurred, that is, 2322 to 2345 
for MIT cell A and 2316 for MIT cell B. For cell A the average area is 14 square boxes, or an 
effective diameter of 3.3 km, close to that noted in section 7.3.2 for 40 dBZ. For times and 
cells with lower peak intensities, the areas are generally much larger and of little interest. 

From the tabulated areas fitted to each contour of log Z it is easy to obtain the apex or 
log Z at zero area. The peak values found are listed in table 7-2 along with the precise peak 
values from the MIT-LL tape, converted to dBZ, It may be seen that even though the loga- 
rithmic data are quantized in 10-dBZ steps, the use of several levels and a good model has 
more effect on accuracy than the quantizing interval. In other words, the model is a suffi- 
ciently good representation of the cell that a peak at 46 dBZ can be predicted from inputs 
at 10, 20, 30, and 40 dBZ. A detailed plot using 1-dB steps shows that the large difference 
at 2334 is due to the coincidence that a notch in the distribution fell at 40 dBZ, shifting the 
best-fit profile using that point. 
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Table 7-2.— Peak Values of Log Z 


Time 

Derived, dBZ 

Precise, dBZ 

Difference 

2322 

47 

46.4 

0.6 

2328 

45.5 

45.8 

-.3 

2334 

a 47.5 

50.0 

-2.5 

2340 

44 

44.3 

-.3 

2345 

42 

40.9 

-1.1 


a Found using point at log Z = 5. 


7.5 RELATED SETS OF ELLIPSES 


Tables 7-3 to 7-5 give the results of ellipse fitting operations. The following information 
is presented in this section with regard to related sets: 


Data source 

Raw data 

Ellipse data 

Data plots 

MIT logarithmic 

Figure 7-1 

Table 7-3 

Figures 7-4 and 7-5 

MIT linear 

Figure 7-1 

Table 7-4 


NSSL 

Figure 7-2 

Table 7-5 

Figure 7-2 


Figure 5-3 has presented horizontal cross sections through MIT-LL cell A, supporting the 
use of a vertical, cylindrical cell model. Figure 7-6 shows vertical sections through MIT-LL 
cells A, B, and C along lines shown in figure 3-4(c). The tapered tops and split bottoms are 
evident in cells B and C. The distributions of the parameters of the ellipses, TAV and TRZ, 
are shown in figures 7-3(c) and 7-3(d). The histogram for NSSL ellipse set orientation, TAV, 
shows the predominant northeast-to-southwest direction assumed by Rogers and Rao (ref. 1 ), 
but the MIT-LL results do not. However, the latter storm cells are weaker, with only a single 
box with log Z = 5. They are also nearly round, with an average minor/major axis ratio of 
about 0.75 and about the same value for centroid line slope PRJ. The NSSL storm, on the 
other hand, showed a maximum intensity of log Z = 5 from 1735 to 1820 and was highly 
elliptical, with a larger PRJ. The histogram for centroid line azimuth, TRZ minus TAV, 
shows that the line of ellipse centers was displaced about 20° ccw from the major axis line in 
data from all sources. 

The NSSL time sequence of figure 7-2 and additional data from table 5-4 are summa- 
rized in figure 7-7. The azimuth (AZ) and range (RU) coordinates can be seen to vary 
smoothly, with the azimuth fairly constant and the range increasing slowly in a southeast- 
erly direction at about 12 km s -1 (23 knots). The average orientation, TAV, runs northeast 
to southwest, as assumed by Rogers and Rao (ref. 1 ), except for one point. The angle be- 
tween the centroid line, TRZ, and TAV remained small until the cell split a second time. The 
area profile intercept dBZ max fell slowly and the slope remained fairly constant, although 
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the area for log Z- 4, AREA (4), dropped sharply after the split. The weakening of the 
storm is also indicated by decreasing ellipticity (increasing RAV) and increasing symmetry 
or concentricity (decreasing PRJ). The intensity, long life, and major axis of about 37 km 
(20 miles) indicate that this may be the type of supercell associated with tornadoes and/or 
hailstorms. 
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Table 7-3.— MIT-LL Logarithmic Ellipse Data, (a) Cell A 
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C F A CFR CF Z TR Z PRJ TAV RAV AZ-CG-RU AREA ER AZ-CG-RU AREA ER AZ-CG-RU AREA ER AZ-CG-RU AREA ER AZ-CG-RU AREA ER 

3C.3 29.5 162 2C 3C.0 29.8 118 21 33.0 31.0 9 71 



Table 7-3.— (c) Concluded 
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r^E THIS yAJ KAT THTA VtJ t * I N RAT THT A RAJ RIN RAT Tt-TA RAJ RIN RAT THT A RAJ RIN RAT 

5357 5C. 7.9 A.l C.52 53. 7.0 3.5 0.59 58. 5.8 2.3 0.39 57. 5.6 l.l 0.25 31. 2.0 0.7 C.36 

cfa cfr cfz tkz prj tav rav az-cg-rl area er az-cg-ru area er az-cg-ru area er az-cg-ru area er az-cg-rg area er 

35. C 35.8 l C 5 3C 35. C 35.8 77 51 35 .5 35.5 53 65 35.1 35.1 17 82 

35.7 35.5 2.5 216 C.7 59 0.3 35.1 35.9 1C7 3C 35.8 35.6 77 53 35.5 35.5 57 69 35.2 35.2 17 86 
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TIPt THT« P A J PIN RAT THTA PAJ PIN RAT TP T A PAJ PIN RAT TPTA PAJ PIN RAT THTA MAJ MIN RAT 
2345 AS. 4.1 3.9 0.95 85. 4.1 3.6 C.86 28. 2.3 1.6 0.68 C. 0.5 0.5 1.C0 0. 0.5 0.5 l.CO 
C F A CFR C F 2 TR2. PRJ TAV RAV AZ-CG-RL AREA ER AZ-CG-RU AREA ER AZ-CC-RIJ AREA ER AZ-CG-RU AREA ER AZ-CG-RU AREA ER 

29.6 41.1 52 40 29.8 41.2 48 42 28.5 39.6 12 73 28. C 42.0 1 0 
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TIRt THTA PAJ MIN RAT THTA MAJ MIN RAT THTA MAJ MIN RAT THTA MAJ MIN RAT THTA MAJ MVN RAT 
1750 5.18.3 3.6 0.19 7.16.8 3.0 0.18 9.12.2 Z.3 0.18 13. 7.9 1.7 0 .22 0. 1.1 0.5 O.AA 
CEA CFK CHZ TRZ PR J TAV RAV AZ-CG-RU AREA ER AZ-CG-RU AREA ER AZ-CG-RU AREA ER AZ-CG-RU AREA ER AZ-CG-RU AREA ER 

59. A 27.5 208 15 60.0 27.1 162 18 62.3 27.0 90 23 6A. 3 27.0 A A A8 
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Table 7-5.— Concluded 
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Figure 7-2.— NSSL data and plots, (a) Time 1735. (b) Time 1750. 
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Figure 7-2 (continued).- (c) Time 1805. (d) Time 1820. 
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Figure 7-2 (concluded).— (e) Time 1835. 
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Figure 7,3.— Histograms of error measures, (a) Individual errors, (b) Set minus individual errors, (c) TAV. 

(d) TRZ minus TAV. 
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Figure 7-4.— MIT-LL cell A plotted individual-ellipse data. 
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Figure 7-5.-MIT-LL cell A plotted ellipse-set data. 
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Figure 7-6.— Vertical sections through MIT-LL cells A, B, and C. 
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Figure 7-7.— Temporal characteristics of NSSL cell B. 
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Section 8 
CONCLUSIONS 

It has been shown that intense storm cells can be usefully represented by digitized 
weather data with a resolution of 1 to 2 km. Any deficiency in such resolution is partly 
compensated by the interpolation and extrapolation from adjacent data using the models 
summarized below. This is important in reducing the quantity of input data required and 
the number of parameters needed by the model. Height slices 2 km apart and a repetition 
rate of four scans per hour seem adequate, but this judgment may be based on inadequate 
sampling. Linear stratification in Z instead of log Z is useful for detailed storm core studies 
but requires a large number of intensity levels in the input data unless only the most intense 
cells are of interest. This may not be the case as higher frequencies come into use, which are 
subject to greater attenuation by the lower intensity areas within a storm. 

Using the above data and contours of constant reflectivity, it was found that— 

( 1 ) Individual horizontal contours of log Z are fairly well fitted by ellipses. 

(2) Sets of contours for different reflectivities at a given time and height are very well 
fitted by assumption of constant ellipse orientation, ellipticity, and uniform spac- 
ing of ellipse centers on a straight line, about 20° ccw from the major axes. 

(3) The ellipse areas in a set were found to vary directly with log Z, in agreement with 
results of other observers. 

(4) For related sets at different heights or lateral positions, fitted ellipse parameters 
were found to vary slowly and uniformly except for the effect of vertical wind 
shear, so that the best model appears to be nested, vertical elliptical cylinders. 

Too few cells were studied to generalize the results, but an effective core diameter of about 
3.5 km was indicated with a favored ellipse orientation of northeast to southwest for the 
strongest storm studied. 

A comparison was made of the attenuations along paths through the centers of five 
horizontal sections of MIT-LL cell A as observed and as modeled by ellipse sets. The agree- 
ment at three of the levels was excellent, but at two levels it was degraded because the ori- 
entation of the core log Z = 4 ellipse remained relatively constant with height, whereas the 
average rotated ccw. It may therefore improve the model for some purposes to use the core 
orientation for the ellipse set instead of the average of the individual orientations. 

Many more storm cells should be processed to verify the model parameters derived and 
to suggest appropriate methods for treating splitting cells, or at least to indicate the magni- 
tude of the problem. 
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